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ABSTRACT 


An aerosol optical depth retrieval algorithm in the visible wavelengths for the 
NOAA POES AVHRR and GOES-8 visible imager is presented for the cloud free, 
marine atmosphere. The algorithm combines linearized single-scatter theory with an 
estimate of surface reflectance. Phase functions are parameterized using an aerosol size 
distribution model and the ratio of radiance values measured in channels 1 and 2 of the 
AVHRR. Retrieved satellite aerosol optical depth (AOD) is compared to three land- 
based sun photometer stations located on islands in the western Atlantic during July and 
September, 2001. GOES-8 channel 1 (visible wavelength) radiance values were initially 
calibrated using techniques developed by Rao. Additional corrections to the channel 1 
GOES-8 radiances were made by applying a linear offset factor obtained during the 
experimental time period through comparison with AVHRR radiances. The results for 
the GOES-derived AOD compare favorably to the AERONET-measured AOD values. 
Eor both NOAA and GOES data, the comparison dataset has a correlation coefficient of 
0.67 with a standard error of 0.07. Eor higher AOD cases (d = 0.25), the general trend 
was for the satellite-derived AOD values to underestimate AERONET-observed 
conditions. During these higher conditions, the scattering phase function pattern 
contained within the algorithm deviated from the expected pattern, especially between 
140° - 180°. Overall, the more accurate calculations of AOD occurred over scatter angles 
between 140° - 150° and 170° - 180°. 
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1. INTRODUCTION 


Monitoring tropospheric aerosols on a global scale is essential for evaluating the 
earth’s radiation budget. Aerosols are known to cause a net cooling effect by scattering 
incident solar radiation back to space and by interacting with clouds in a way that 
increases overall albedo (Charlson et al., 1992, Twomey et al., 1977, IPCC, 1996). King 
et al. (1999) reports that the impact of aersosol radiative forcing, both directly 
(scattering) and indirectly (interaction with clouds), produces a cooling range of 
-1.4+1.5 W rri^ to -2.5+2 W m‘^. This result offsets the well-known concept of the 
greenhouse warming impact, estimated to be +2.5+0.3 W m' . 

Assessing aerosol properties also has military implications. Infrared wavelength 
image and ranging systems, laser-guided weapons systems, and laser communication 
systems are sensitive to the environment. Infrared imaging and ranging systems are 
strongly affected by variations of atmospheric aerosols. Laser systems used for radio and 
satellite communications operate within the visible and near-IR wavelength ranges, and 
can be greatly affected by varying aerosol properties, such as density, size distribution, 
chemical and physical composition (Bloembergen et al. 1987; Cordray et al. 1977). 
Brown (1997) reports that the proper interpretation of aerosol radiative properties in the 
coastal zone is “important to the design, planning, and operation of electro-optical 
weapons and sensor systems near coastal boundaries”. 

Another research area of global aerosol impact focuses on the transport of dust and 
pollutants from one region to another. ACE-Asia is a 4 year project (2001 - 2004) 
devoted to the study of aerosol profiles in the Pacific basin generated by desert dust and 
industrial pollution over Asia. The Puerto Rico Dust Experiment (PRIDE) in 2000 
studied the impacts of African desert dust that is transported over t) the Caribbean and 
the eastern U.S. 

Given the challenges listed above, there is a developing interest to globally quantify 
aerosol properties on fine spatial and temporal scales. Thus far, this analysis has proven 
to be a daunting task, since most established aerosol sensing is land-based, providing 
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poor spatial and temporal coverage. Higurashi et al. (1999) suggests that aerosol 
concentration, size distribution, composition, and optical properties will have to be 
measured globally, and that satellite remote sensing is an effective tool for such a task. 
Over the past few decades, scientists have developed algorithms to convert satellite 
upwelling radiances into aerosol properties such as optical depth. So far, most of the 
aerosol remote sensing studies have used the NOAA Advanced Very High Resolution 
Radiometer (AVHRR) channel 1 and 2 sensors. These algorithms were developed by 
assuming certain aerosol characteristics before processing the upwelling radiances 
(Durkee et al. (1991), Kaufman et al. (1990), Higurashi and Nakajima (1999)). 

This study focuses on one such algorithm developed by Durkee et al. (1992) 
(hereafter referred to as the NFS algorithm). The NFS algorithm ingests AVHRR data 
within a cloud-free, single scatter environment. By using the ratio of channel 1 and 2 
radiances, an estimate of the aerosol size distribution is extracted. During three recent 
field campaigns, Durkee et al. (1999) has shown that the NFS algorithm performs well 
for aerosol optical depth (AOD) below about 0.4 at 0.63|im wavelength. But the results 
only provided snapshots of the experimental regions, since AVHRR passes over a 
particular region a few times per day. Brown (1997) incorporated both the AVHRR and 
the Geosynchronous Operational Environmental Satellite (GOES) satellite to provide 
temporal coverage of AOD over an experimental region. In the upgraded NFS algorithm, 
an aerosol model index is derived from the AVHRR data is applied to the GOES 
retrieval. The results were encouraging, but limited. The focus of this report is to expand 
this approach. The following issues will be addressed: 

• Froper radiance calibration of the visible sensor of GOES 

• Validation of the AOD derived from NOAA and GOES-8 

• Evaluation of the phase function parameters used in the retrieval algorithm 

Chapter II describes the radiative transfer theory and the simplifying atmospheric 

assumptions used in the satellite optical depth retrievals. Chapter III describes the data 
sets and the instrumentation used. Chapter IV describes the calibration and correction 
techniques applied to the retrieved GOES channel 1 radiance. In addition, the AOD 
retrieval procedures for both the AVHRR and GOES are discussed. Chapter V discusses 
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the results and Chapter VI presents the conclusions and recommendations for future 
research. 
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11. THEORY 


A. RADIATIVE TRANSFER THEORY 

Estimates of AOD from satellite-sensed upwelling radiances are calculated by the 
Radiative Transfer Equation (referred to hereafter as RTE). The RTE takes into account 
the various paths and interactions that solar radiation encounters with atmospheric 
particles and the earth’s surface before reaching the satellite sensor. Erom Liou (1980), 
the general form of the RTE for a given solar wavebngth in a plane parallel atmosphere 
is: 


9 


dL(d,0) 

dd 



(o,o'Jdo' 


( 1 ) 


where: 


L = diffuse radiance (WW pm sr) 
5 = optical depth, 


(Do = single scatter albedo, 
p = COS0 (0 = satellite zenith angle), 

Po = COS0O (00 = solar zenith angle) 

Q = solid angle (0,(p) (cp = azimuth angle), 
p(Q,Qo) = scattering phase function and 
tiEo = incoming solar radiative flux. 

The terms on the right hand side of the equation describe the radiation loss due to 
atmospheric attenuation, the intensity gained by multiple scattering into the beam and the 
beam addition due to single scatter events. 
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For the aerosol optical depth environment, Brown (1997) and Durkee et al. (1991) 
greatly simplify the Eq. (1) by making the following assumptions: the satellite scans the 
region only in a cloud-free, marine environment, with a single scatter approximation. In 
addition, contributions from ozone absorption, Rayleigh and ocean surface effects are 
eliminated from the calculation. The equation is then reduced to: 


L 

a 



( 2 ) 


where the subscript “a” refers to aerosol related quantities and 'Fs represents the 
scattering angle. In Eq. (2), the incoming solar radiance (Eo) is a known constant and the 
satellite zenith angle (p) is obtained by the satellite-earth geometry. The single scatter 
albedo is a measure of the ratio of radiance scattered versus radiance absorbed. Eor 
marine aerosols sensed by visible wavelengths, the particles (salt, sulfate) are weakly 
absorbing and therefore, the single scatter albedo (cOo) is nearly one. La represents the 
satellite-detected upwelling radiance only from contributions of atmospheric aerosol 
particles. Radiances due to sun glint (Cox and Munk, 1954), Rayleigh scatter and surface 
reflected radiance were eliminated. Brown (1997) describes in detail the steps involved 
to eliminate these radiances. This leaves the single scatter phase function, p('Fs), and the 
AOD, 5a, which are discussed below in greater detail. 

B. OPTICAL DEPTH 

Solar radiation interacts with the atmosphere in three ways: absorption, scatter 
and emission. The density and absorbing properties of the particles determine how much 
incoming solar radiation is attenuated before reaching the satellite sensor. The 
combination of absorption and scattering is called extinction, which is described by the 
following equation: 


Sext = jpr^Qext (m,r)n(r)dr (3) 

0 

where Oext is the extinction coefficient, which describes how easily the particle interacts 
with the solar photons. Eor the remaining terms, r is the particle radius, Ttr^ is the particle 
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cross-section area, Qext(ni,r) is the extinction efficiency factor that varies with the radius 
(r), wavelength and composition, m is the complex index of refraction, and n(r) is the 
number of particles in a size increment r to r -i- dr. Changes in the size, composition, and 
distribution of constituents of suspended particles in the atmosphere lead to a change in 
extinction. Therefore, measurement of extinction can provide information about the 
characteristics of the atmosphere’s particulate distribution. The scattering of solar 
radiation is the largest contributing factor to extinction in the visible and near infrared 
wavelengths. Since this report deals with aerosol extinction only, Rayleigh scattering can 
be accounted for and factored out of the optical depth calculation because the 
atmospheric molecular particles are nearly uniform both spatially and temporally (Durkee 
et al, 1991). There is also a small contribution of absorption by ozone gases within the 
upper atmosphere, but its effect is minimal and easily factored out as well. This leaves 
scattering effects due only to aerosol particles. Scattering due to cloud particles are not 
considered. Therefore, for a cloud-free marine environment, the extinction coefficient 
can be approximated by the scattering coefficient: 

oo 

'‘scat =lP''^‘3scat 

0 

where Oscat is the extinction coefficient due to scattering, and Qscat(ni,X,r) is the scattering 
efficiency factor, representing the ratio of total energy scattered in all directions to 
incident energy. 

Optical depth, as sensed by satellite is expressed as the sum of the atmospheric 
extinction integrated vertically from the surface through the atmosphere and is expressed 
as: 

H H 

^ ~ 1 ^ ext^^ ~ 1 ^ scat^^ 

0 0 

where “H” represents the height of the satellite above the earth’s surface. As mentioned 
earlier, the study is conducted under atmospheric conditions where most of the extinction 
is due to scatter. 
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C. SCATTERING PHASE FUNCTION (p(\l/s)) 

The scattering phase function determines which direction the radiation is scattered 
when it encounters a particle. In the case of atmospheric aerosol particles, the satellite 
senses solar radiation that encounters a particle and is scattered into the viewing direction 
of the sensor. This scatter angle is illustrated in Figure 1. The scattering phase function 
is dependent on the size, shape and distribution of the atmospheric particles. In the case 
of aerosol particles, Mie theory is used to quantify the phase function. Knowledge of 
aerosol size distribution and composition is a necessary step in the proper assessment of 
the phase function. However, identifying the actual aerosol size distribution and 
composition over a particular region is nearly impossible. Therefore, the phase function 
must be parameterized. Durkee et al. (1991), uses the ratio of the two visible wavelength 
channels of the NO A A AVHRR since the ratio is sensitive to changes within the aerosol 
size distribution. 


Using a single scatter model, there are three possible paths by which solar 
radiance can interact with aerosol particles and eventually reach the satellite sensor. 
Figure 2 displays these paths. The first path (path (a)), known as direct scatter, occurs 
when incident solar radiation scatters off of an aerosol particle, in a backscatter direction, 
directly to the satellite sensor. The other two paths are created by reflectance off of the 
ocean surface. Path (b) is solar radiation reflected off the ocean surface and then forward 
scattered (0 - 90°) off of an aerosol particle and toward the satellite sensor. Path (c) is 
forward scatter (0 to 90°) off of an aerosol particle, and then is reflected off the ocean 
surface and toward the satellite sensor. The reflectance off of the ocean surface can be 
estimated using Fresnel reflection coefficients, as shown in the following equation: 


U = 0.5 


V 


sin(./i- 


where: 


0i: angle of reflectance 


( 6 ) 
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0t: angle of transmission 

From Snells’s Law, 0t = sin'^(sm(0i)/m)) where m is the index of refraction of seawater 
(1.33). Combining all three scatter paths, the result is an effective phase function that is 
defined by: 

Peff =P-+P+(rM„+rM) (7) 

where: 

Peff = effective phase function 

P- = phase function at back scatter angle (90 to 180°) 

P+ = phase function at forward scatter angle (0 to 90°) 
r = Fresnel reflection coefficients at both po and p 

It is the effective phase function value (peff) that is applied to the optical depth calculatbn 
in Eq. (2). 

D. SUMMARY OF ASSUMPTIONS TO THE NPS ALGORITHM 

For this study, the NPS algorithm processes AOD output by applying a simplified 
version of the radiative transfer theory, shown in Eq. (2). This process was accomplished 
by noting the following assumptions. The algorithm was applied to regions free of cloud 
and sun-glint. The aerosols are assumed to be spherical and non-absorbing, and the 
optical depth is low enough to be dominated by single scattering. These assumptions 
have been shown to provide reasonable results within low to moderate AOD conditions 
(Brown, 1997 and Smith, 1998). 
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Figure 1. 



Polar plot of scattering phase function describing the scatter angle of incident 
radiation with an aerosol particle. 



Figure 2. Schematic of solar radiation trajectories that interact once (single scatter) with 
aerosol particles and eventually reach the satellite sensor. Path “A” describes direct 
scatter while paths “B” and “C” indicate diffuse scatter. 
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m. DATA 


To validate the optical depth retrieval method, case study days were chosen based 
on the availability of matching datasets between GOES, NOAA and AERONET. This 
chapter briefly describes the data sets and instrumentation used for this study. 

A. INSTRUMENTS 

1. NOAA Advanced Very High Resolution Radiometer (AVHRR) 

The AVHRR instrument senses upwelling radiances of 5 channels, ranging from 
visible to infrared. Table 1 describes the bandwidths of these channels. The AVHRR 
instrument is part of the National Oceanic and Atmospheric Administration (NOAA) 
Polar Orbiting Operational Environmental Satellite (POES) series of satellites. These 
satellites are in sun synchronous orbit at a height of 883 km and provide at least four 
passes per day over a given region of the earth. The AVHRR scans at nadir with a width 
of approximately 2000 km and a sub-satellite pixel resolution of 1.1 km by 1.1 km. 

The AVHRR instruments onboard the NOAA-14 (launched on 01 January 1995) 
and NOAA-16 (launched on 02 Eebruary 2001) satellites provided the data for this study. 
Over the experimental region, the NOAA-16 AVHRR provided the local afternoon data 
while the NOAA-14 AVHRR provided data late in the afternoon. The data from the 
NOAA-14 was at times questionable due to low sun angle problems. Therefore, the 
NOAA-16 data was the more reliable dataset. 

The AVHRR dataset was transmitted to the receiver at Wallops Island, Virginia, 
and archived by the NOAA National Environmental Satellite, Data and Information 
Service (NESDIS) Satellite Active Archive, in Eevel lb format. AVHRR Eevel lb data 
is in 10 bit precision format that have been quality controlled, assembled into discrete 
data sets, and to which earth location and calibration information has been appended, but 
not applied to the data. Other parameters included are time, quality flags, solar zenith 
angles, and telemetry. All AVHRR channels are calibrated prior to launch. Channels 1 
and 2 have no onboard calibration systems. Post calibration methods for these channels 
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were devebped by the NOAA/NESDIS Office of Research Applications based on results 
of Rao and Chen (1995).CHECK SPACING HERE, BETWEEN TEXT AND HEADING 

2. GOES-8 Imager 

The GOES-8 imager is part of the Geosynchronous Operational Environmental 
Satellite (GOES) series of satellites that covers the western Atlantic region. GOES-8, 
launched on 13 April 1994, operates at a geosynchronous orbit at 36,000 km over the 
equatorial subpoint at 75 W. Eike the AVHRR, the GOES imager scans in 5 channels. 
However, the GOES only senses in one visible wavelength. In addition, the scanning 
resolution varies between channels. Table 2 provides the bandwidths of the GOES 
imager. Eor this research, only the visible channel (channel 1) data is used for the AOD 
calculations. Eike the NO A A AVHRR, the visible channel does not have an onboard 
calibration system. In addition, immediately after the launch, there was a nonlinear 
decrease in signal strength. Compared to the AVHRR, the GOES dataset was much 
noisier and required several steps of adjustments before the data was comparable to the 
expected radiances. Chapter 4 describes these correction procedures. 

3. AERONET Sun-sky Scanning Spectral Radiometer 

Data collected from the Aerosol Robotic NETwork (AERONET) automated 
radiometers are applied as ground truth of AOD for this study. AERONET contains a 
global network of ground-based automated radiometers, each of which consists of a 
CIMEE Electroniques 318A spectral radiometer that makes direct sun or sky 
measurements. The direct sun measurements are made in 8 spectral bands; 340, 380, 
440, 500, 670, 870, 940, and 1020 nm. The 670 nm wavelength data is comparable to 
that used in the optical depth calculation obtained from both the GOES and NOAA 
retrievals. Optical depth is calculated from the spectral estimation of direct beam 
radiation based on the Beer-Bougner Eaw. As with the NPS AOD algorithm, attenuation 
due to Rayleigh scatter and absorption of ozone is taken into account. In addition, 
contributions from gaseous pollutants are also eliminated. Detailed information about the 
operations and instruments of the AERONET system is provided in Holben et al. (1998). 

AERONET data was obtained from radiometer instruments installed on three 
island sites: eastern Bermuda, (U.K.) (32 22’N/64 41’W), Ea Paguera, Puerto Rico (17 
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58’N/67 02’W) and Guadaloup, Island (Fr.) (16 19’N/61 30’W), all within the western 
Atlantic Basin (see Fig. 3). For all 3 sites, the dataset contains only the times that have 
not been contaminated by clouds. Data from the La Paguera site has an additional quality 
control check for clouds. 


Table 1. NOAA AVHRR Radiometric Channels. 


Channel 

Band Widths (|ini) 

1 (visible) 

0.58-0.68 

2 (Near Infrared) 

0.725- 1.10 

3 (infrared) 

3.55-3.93 

4 (infrared) 

10.3-11.3 

5 (infrared) 

11.5-12.5 


Table 2. GOES Imager Radiometric Channels 


Channel 

Band widths (|im) 

Resolution (km) 

1 (Visible) 

0.55-0.75 

1 

2 (infrared) 

3.80-4.00 

4 

3 (infrared) 

6.50-7.00 

8 

4 (infrared) 

10.20-11.20 

4 

5 (infrared) 

11.50-12.50 

4 
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IV. METHODOLOGY 


A. SELECTION OF CASES 

Datasets were collected on a daily basis from the GOES-8, NOAA, and 
AERONET instruments. The process started in the middle of July 2001 and lasted until 
the end of September 2001. GOES-8 data was collected from a Terascan receiving 
platform located at the Naval Research Eaboratory, in Monterey, CA. The data was 
downloaded in a short (10 bit) format to a UNIX workstation every half-hour. GOES-8 
data in August was downloaded into 8 bit instead of the 10 bit format, and was not 
suitable for processing. As a result, all data in August was rejected. 

NOAA data was readily available through a NOAA NESDIS archives. Of the 
two satellite sensors, data from NOAA-16 AVHRR was more suitable than the NOAA- 
14 because NOAA-16 orbital passes occurred over the study region (eastern Atlantic and 
Caribbean Sea) during the times of interest - late morning through early afternoon. 
NOAA-14 orbital passes occurred during the late afternoon and its data was often 
rejected because of low sun angles. Eor both NOAA sensors, datasets were also rejected 
when the sun glint pattern (specular solar reflection) covered the study area. The NFS 
algorithm identifies and rejects the contaminated region during data processing. 

AERONET data was also readily available via NASA’s AERONET archive for 
all 3 locations. Cases were rejected if there was insufficient data covering the particular 
site during the local afternoon hours. These conditions usually occurred because of the 
presence of clouds over the AERONET radiometer’s field of view, which unfortunately, 
is a common diurnal occurrence over all 3 islands, especially during the afternoon hours. 
Data was manually rejected during times that the observed AOD values showed 
tendencies of significant increase before a cloudy episode. 

B. GOES-8 CALIBRATION AND CORRECTION PROCESSES 

Although extensive studies have shown that the NFS algorithm performed well 
with NOAA AVHRR data, the results using GOES-8 data were very limited. Eor GOES- 
8, there is no on-board calibration for the channel 1 radiances. In addition, there was a 
problem of non-linear signal degradation immediately after the launch in 1994. There 
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have been several attempts to perform vicarious calibration techniques to adjust for the 
weakening signal strength and to take into account the post-launch degradation. There is 
an additional problem with noise and signal radiometric resolution generated from GOES 
sensors. Figure 4 displays an example of radiance comparisons between the calibrated 
GOES and AVHRR datasets. The image passes occurred at a similar time and were 
registered over the same 100 km by 100 km domain surrounding Bermuda. In addition, 
both satellites had similar viewing geometries (scatter angles were + 0.1° of each other). 
The atmospheric conditions immediately south of Bermuda (outlined in red) were clear 
and homogeneous at this time. The panel to the right of Figure 4 are plots of the 
histogram frequency distribution of radiances that were extracted within the outlined 
region. As expected, the NOAA-16 image in (b) displays a relatively homogeneous field 
of radiances; its corresponding frequency histogram displays shows a pronounced signal 
peak with a narrow radiometric width, indicative of the pristine atmospheric conditions. 
In contrast, the GOES-8 sensor, situated in an orbit that is 40 times the distance of the 
NOAA sensor, produces an image (shown in (a)) that is significantly noisier; its 
corresponding histogram profile displays a weaker signal peak and wider radiometric 
range. As atmospheric conditions become hazier, the GOES-8 peak signal and 
radiometric resolutions become even less discernable, thus complicating the processed 
AOD calculations. In addition, the GOES-8 signal peak at ~ 16 Wm“^sr‘Vm^ is 
significantly weaker than the NOAA-16 signal peak at -22.5 Wm" sr' pm' , thus 
necessitating a further correction factor to GOES-8 for further processing. 

For this study, in order to match GOES with NOAA data during AOD processing, 
two correction techniques were applied to the GOES channel 1 radiance data. Dr. C. R. 
N. Rao (personal communication in July, 2000) developed a calibration methodology of 
GOES-8 channel radiance by a vicarious technique, selecting a radiometrically stable 
calibration site located in the Sonoran desert (34.0°N/114.1°W). Radiometrically stable 
calibration is defined as the long term mean value at the top of the atmosphere albedo that 
remains uniform in time. Details of this method can be found in Rao and Zhang (1999) 
and Rao et al. (1999). 
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Dr. Rao developed a simplified version of the correction factor (gain factor) that 
can be expressed in the following equation: 


GAIN_ FACTOR = = 1.192*(l+1.688E-04*d) (8) 

where: 

d = number of days since the launch date of the GOES-8 satellite 

A(pre) = albedo (%) calculated using the pre-launch calibration coefficients 

A(d;post) = albedo (%) calculated on day ‘d’ after the launch of GOES-8, that 
accounts for degradation in orbit) 

Table 1 in Appendix A lists the calibration factor applied to GOES-8 for each 
case study day in this report. However, a preliminary assessment of the calibrated GOES 
data indicated that its resulting AOD values were significantly higher than the NOAA- 
generated AOD as well as the AERONET observations of AOD. Therefore, a further 
correction method was applied, as discussed below. 

The correction technique involved comparisons between GOES-8 and NOAA-16 
channel 1 radiances, whose wavelengths, centered on 0.65 and 0.63 pm, respectively, 
were similar. The assumption is that the calibrated AVHRR channel 1 radiances are 
accurate. The process involved the collection of a sample of cases when the GOES and 
NOAA viewing geometries over a selected location were similar. GOES and NOAA 
viewing geometries were defined as ‘similar’ when the scatter and azimuth angles were 
within + 0.1° of each other. Eigure 5 presents the comparisons over 296 pairs of GOES 
and NOAA channel 1 radiances. As shown, the radiance ranged along the low end of the 
radiance spectrum (0 to 40 Wm" sr' pm' ) which is where the detection of aerosols would 
occur. As shown, there was very poor correlation between the NOAA and GOES data, 
due to the large noise problem in GOES. As a plausible correction, it was decided to 
perform a manually-determined selection of the “center of mass” within the domain 
shown in Eigure 5. Using the Cartesian coordinates, the selected center of mass of the 
distribution was positioned at point Y (NOAA axis) = 24.0 Wm' sr' pm' and point X 
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-9 -1 -1 

(GOES axis) = 19.0 Wm sr' |am‘ . Assuming there is a linear relationship between the 
channel 1 radiances of GOES and NOAA, the correction factor was determined by 
locating the X (GOES) intercept from the slope (red dashed line), which was estimated to 
be 5.5 Wm' sr" pm' . Therefore, before the actual AOD processing took place, the 
correction ion factor was added to the calibrated value of the GOES channel 1 radiances. 

C. AOD RETRIEVAL FROM SATELLITE DATA 

Eigure 6 displays a flowchart to the steps involved in the processing of fie 
satellite data. A software algorithm (hereafter referred to as the “NFS algorithm”) was 
developed by Brown (1997), and then modified by both Smith (1998). The NFS 
algorithm consists of 3 major parts: pre-processing, processing, and post-processing. For 
the pre and post-process stages, a combination of Terascan, Cshell and FERE software 
manipulate the data. During the processing stage, the satellite data is processed and 
converted to AOD values by code written in FORTRAN 77. 

1. Pre-Processing Stage 

For the preprocessing portion, the NOAA AVHRR raw data is initially converted 
from Eevel lb to Terascan data Format (TDF). The data is then calibrated from raw 
radiance counts into albedo (%) from the visible wavelengths (channels 1 and 2), and 
temperature (“K) within the near IR and IR wavelengths (channels 3, 4, and 5). The 
GOES raw data is already in this format. Both the GOES and NOAA data then have the 
earth-sun-satellite geometry angles attached to the data, which consists of the solar and 
satellite zenith angles, the azimuth angle and the scattering angles. These angles were 
obtained from the telemetry dataset stored within Terascan. Next, the data is registered in 
Mercator coordinates to a predetermined area. For this study, the areas consist of 200 km 
X 200 km domains surrounding Bermuda, Fa Faguera, Fuerto Rico, or Guadaloupe 
Island. The final stage in the preprocessing is the conversion of the data from TDF to 
binary format for the next stage. 

2. Processing Stage 

During the process stage, the radiances are converted to AOD for both GOES and 
NOAA data. For GOES, the channel 1 radiances data is initially calibrated and a 
correction factor is applied following the process described in the previous section. For 
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both datasets, sun glint contamination is removed by applying a method used by Cox and 
Munk (1954). 

AOD is calculated following Eq. (2), 


^a 


(4jiLa)/ 

/(?oFoP(?s))' 


The cosine of the satellite zenith angle (p), the single-scatter albedo (? o), and the solar 
radiance (Fo) are all constants. Radiance due to aerosol scatter (La), is mathematically 
straightforward and is described in detail by Brown (1997). The scattering phase 
function, p(? s), is obtained from a process described below. 

a. Scattering Phase Function Processing 

Obtaining the scattering phase function values requires knowledge of the 
aerosol characteristics and size distribution, which is not routinely available. Therefore, 
the scattering phase function must be parameterized. Durkee et al. (1991) developed the 
parameterization technique used within the NFS algorithm. The technique consists of 
calculating the ratio of the NOAA channel 1 and 2 radiances, ‘S 12 ’. The scattering 
efficiency (Qscat) of an aerosol distribution is wavelength dependent and peaks when the 
radius of the aerosol particle is nearly equal to that of the radiation wavelength. As a 
result, S 12 will be larger for smaller size particle distributions and smaller for larger size 
aerosol particle distributions. Si 2 varies from pixel to pixel. Therefore, variations in the 
aerosol size distribution can be detected within the pixel resolutions of the satellite image 
data. 


Brown (1997) generated seven models of aerosol size distributions. These 
models typify general conditions within the maritime environment. The scattering phase 
function and extinctions for these models were calculated using Mie theory. Table 3 
describes the attributes for each of the 7 models. These distributions consist of one 
single-mode and 6 two-mode log normal distributions with varying radii and standard 
deviations used to describe variations of aerosol distribution widths in the maritime 
atmosphere. Figures 7, 8, and 9 illustrate the effect of the aerosol size distribution 
models on S 12 and the scattering phase function (developed by Brown, 1997). Figure 10 
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is a composite chart that illustrates the actual phase function extraction process. This 
process only applies to NOAA processing. For each pixel, the combination of the scatter 
angle and the computed Si 2 value is entered into a lookup table (LUT) represented in the 
upper left portion of Figure 10 to determine the model aerosol distribution that is 
consistent with the measured Si 2 . An aerosol model index (AMI) is interpolated between 
models MO and M6. In the example shown, the interpolated value is situated between 
Models M2 and M3 (i.e., AMI is approximately 2.5). The AMI values are collected for 
each pixel and stored in a file. During the GOES processing, this data from the AMI file 
is accessed. For both GOES and NOAA processes, the next and final step to the NOAA 
phase function processing is entering the AMI value with the scatter angle to the 
scattering phase function EUT, as graphically represented in Eigure 4.4 and in the lower 
left portion of Eigure 10. Eor GOES processing, the AMI values are assumed to be 
constant during the entire time period. At this stage, all of the input parameters have 
been determined and the AOD is calculated using Eq. (2). 

3. Post-Processing Stage 

During the post-processing stage, all of the original channel data, calculated data, 
and extracted scattering phase function values are reformatted back to TDE; this data is in 
image form and can then be viewed and analyzed via the Terascan visualization software. 
Table 4 lists the output products generated by the NPS algorithm. Eor this study, only the 
channel 1 AOD is calculated for the case study analysis. 
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Figure 4. Comparisons of radiance images and associated frequency of radiance histograms 
between (a) GOES-8 and (b) NOAA-16 on 02 September 2001. Histograms were 
developed from areas within red annotations. Radiances, as shown within the color 
legend and the histogram x axis, are in units of Wm'^scVirf'- 
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Figure 5. Comparisons of Channel 1 radiances between GOES-8 and NOAA-16. 

-9 -1 -1 

Radiances are in units of Wm" sr" pm' . 
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Figure 6. Satellite AOD retrieval process (portions obtained from Brown, 1997). 
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Figure 7. Model aerosol size distributions (from Brown, 1997). 
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Figure 9. Model S12 values (from Brown, 1997). 
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Figure 10. Parameterization of the scattering phase function, p(? s), is described, (a) Aerosol 
model size distributions, (b) size index, S 12 , (c) scattering phase function calculated from 
the model size distributions as a function of scattering angle, (portions of figure obtained 

from Durkee et al., 1999.) 
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Table 3. Characteristics for each of the 7 models of the aerosol size distribution. Mode 1 


models the background aerosols while Mode 2 models the ocean-produced aerosols. 


Model 

Mode Radii (pm) 

Mode 1 Mode 2 

Number density (N) 

Mode 1 Mode 2 

Std dev (s) 

Mode 1 Mode 2 

MO 

0.1 

0.0 

1000 

0 

1.7 

0.00 

Ml 

0.1 

0.3 

1000 

3 

1.7 

2.10 

M2 

0.1 

0.3 


5 

1.7 

2.20 

M3 

0.1 

0.3 

1000 

8 

1.7 

2.35 

M4 

0.1 

0.3 

1000 

10 

1.7 

2.51 

M5 

0.1 

0.3 

1000 

13 

1.7 

2.60 

M6 

0.1 

0.3 

1000 

15 

1.7 

2.70 


Table 4. List of the output products for each pixel within the image file as generated by the 


NFS algorithm for both GOES and NO A A data processes. 


GOES-8 products 

NOAA-14 and 16 products 

Channel 1 total radiance 

Channel 1 total radiance 

Channel 2 total radiances 

Channel 1 Rayleigh radiance 

Channel 1 Rayleigh 
radiance 

Channel 2 Rayleigh radiance 

Channel 1 aerosol radiance 

Channel 1 aerosol radiance 

Channel 2 aerosol radiance 

Satellite zenith angles 

Satellite zenith angles 

Solar zenith angles 

Solar zenith angles 

Relative azimuth angles 

Relative azimuth angles 

Scatter angles 

Scatter angles 

Channel 1 phase functions 

Channel 1 phase functions 

Channel 2 phase functions 

Channel 1 AOD 

Channel 1 AOD 

Channel 2 AOD 


Ratio of channel 1 and channel2 (SI2) 


Aerosol Model Index (AMI) 
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V. RESULTS 


For this study, 22 cases were selected to analyze the performance of the NFS 
algorithm in providing sate Hite-derived AOD calculations. In this chapter, Section 5.A 
discusses the results in a clear environment and Section 5.B. discusses conditions within 
a dust environment. Finally, Section 5.C summarizes all 22 cases. 

A. CASE 25 SEPTEMBER 2001. LOW AOD CONDITIONS OVER 

BERMUDA 

1. Synoptic Discussion 

Figures 11a and b present composites of true color images generated by the polar 
orbiting sun synchronous Sea-viewing Wide Field of view Sensor (SeaWiFs). SeaWiFs 
is operated and managed by NASA’s Goddard Space Flight Center and orbits the earth at 
an altitude of 705 km with a 1.1 km resolution at nadir. Although SeaWiFs is primarily 
designed to sense ocean surface properties, a benefit to meteorological research is its 
ability to detect atmospheric aerosol characteristics at high spatial resolutions. As shown 
in Figures 11a and b, dust is quite visible west of Africa, and extends westward into the 
center of the Atlantic Ocean basin. The region surrounding Bermuda is within a clear 
environment, situated north of the dust and east of sulfate concentrations. 

Figure 12 displays dust and sulfate concentrations as generated by the Navy 
Aerosol Analysis and Prediction System (NAAPS) (see Westphal, 2002) on 25 
September 2001 at 18:00 UTC. In this figure, the NAAPS image was formatted to 
display both dust and sulfate concentrations with minimum optical depth thresholds of 
0.1. As in Figure 11, Figure 12 depicts a clear environment surrounding Bermuda during 
the time of this case study. The 850 mb wind field pattern in Figure 12 suggests that the 
region surrounding Bermuda will remain clear during the immediate time period. 

Figure 13 is a visible image of GOES-8 at 17:15 UTC, which is situated about 
local noon in the center of the image. As shown, the island of Bermuda is situated south 
of an eastward moving intense tropical depression. Bermuda is also located east of a 
strong cold front moving off of the US Atlantic seaboard. The immediate region 
surrounding Bermuda is within light convective activity. The associated scattered clouds 
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remain around the island throughout the day, but did not impact the AOD calculations 
conducted over water just east of the island. 

2. AOD and Phase Function Analysis 

Figures 14 and 15 show a series of AOD images processed by the NFS algorithm 
from GOES-8 data. The images occur every 60 minutes from 13:15 UTC to 20:45 UTC. 
The bright white and adjacent yellow features depict small-scale cumulus clouds 
surrounding Bermuda. The general aerosol conditions around Bermuda have very low 
AOD values, ranging from 0.0 to 0.1. Throughout the time range, clear conditions persist 
just east of the island. The location of the AERONET station is on the eastern edge of 
the island. The annotated boxes represent the locations where histogram measurements 
of center of mass took place. The boxes were manually selected and have approximate 
dimensions of 20 km x 20 km. Throughout the time period, the average distance between 
the AERONET station and the center of the box was approximately 25 km. As shown, 
clouds are always nearby, but do not interfere with AOD measurements within the boxes. 

As discussed in Chapter 4, the GOES-8 data is quite noisy. In contrast to Eigures 
14 and 15, the NO A A-generated AOD image in Eigure 16 is much smoother (less noise). 
Extracting reasonable AOD values from GOES and then comparing these values directly 
to AERONET observations was not a straightforward process. The following approach 
to both GOES and NOAA-generated AOD provided reasonable estimates to the study. 

The method of measuring AOD is presented in an example within the frequency 
histogram plot in Eigure 17. A curve is manually drawn (as shown by a solid red curve) 
representing the the frequency distribution of AOD. Then, the mode of the distribution 
was determined. As shown in the figure, the mode for this example translates to an AOD 
value of 0.03. The values on either side of the AOD value (?di: -0.11 and -1-0.15) 
describe the variability about the AOD estimate. As mentioned earlier, GOES sensing 
generally has a wider range of radiances than NOAA data. These values were obtained 
by drawing a horizontal bar at a level that is half of the maximum peak of the normalized 
curve. Eigures 18, 19, and 20 display all of the corresponding sets of histograms for each 
of the red boxes annotated within the images presented in Eigures 14, 15, and 16. 
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Figure 21 presents the time series plot of AOD obtained from Figures 18, 19, and 
20. There is one NOAA measurement obtained at 18:16 UTC. “Variability” bars are 
plotted for each measurement. Because of the inherent noise within the GOES-8 
radiances, the variability of AOD is significantly higher for GOES-8 (?d: +/- 0.10 to 
0.15) than NOAA data (?d: +/- 0.02). Although AOD values can never fall below zero, 
the first two GOES points (Time = 13:15 UTC and 13:45 UTC) produce slightly negative 
AOD values, due again to the noise associated with GOES data. As shown, there is good 
agreement between the AERONET and GOES-derived AOD after about 17:00 UTC. 
The AVHRR-derived AOD value of 0.12 is significantly greater than the corresponding 
AERONET observed value of -0.05 at 18:16 UTC. Included within Eigure 21 is the 
scatter angle pattern shown in a dashed line. Because of the geometric configuration 
between the sun and GOES-8 satellite, local noon occurs approximately where the scatter 
angle position is at its peak, which is also the closest to a direct backscatter configuration 
between the sun and satellite sensor. As time increases, the GOES-8 - sun geometry 
results in more of a side scatter. This geometry affects the phase function determination, 
which is discussed below. 

Eigure 22 displays the phase function values generated from both the NFS 
algorithm and AERONET observations. The satellite-derived phase function values 
(blue diamonds) are obtained from the algorithm’s lookup table, as described in 
Chapter 4. The corresponding red square at each scatter angle represents the phase 
function value required for the satellite-derived AOD to match AERONET AOD. As 
shown, all values are within the backscatter portion (90° - 180°). This figure provides an 
evaluation tool to determine the proper phase function pattern, given the aerosol 
conditions for a particular case. In this case, the AERONET phase function values are 
consistently lower than the satellite-derived phase functions. 

B. CASE 18 SEPTEMBER 2001. HIGH AOD CONDITIONS OVER 

GUADALOUPE ISLAND 

1. Synoptic Discussion 

During the summer months, dust generated from the African deserts are often 
propagated across the southern latitudes of the Atlantic Ocean basin by the easterly trade 
winds, oftentimes impacting the visibility and aerosol characteristics over regions of the 
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Caribbean and the east coast of the US. For this case, the SeaWiFs images in Figures 23a 
and b display several pockets of dust plumes, one of which is approaching the island of 
Guadaloupe on 17 September (Figure 23a). One can surmise that the plume has crossed 
over Guadaloupe a day later on 18 September (Figure 23b). On 18 September, the 
scattered small cumulus clouds surround the Guadaloupe Island region. 

The NAAPS model in Figure 24 shows the leading edge of a significant dust 
plume (AOD > 1.0) reaching Guadaloupe on 18 September at 18:00 UTC. Contours of 
very high optical depth values (AOD > 0.4) are shown within this plume. ^ The 850 mb 
wind pattern at 18:00 UTC supports the continuing propagation of the dust plume over 
the Guadaloupe and Caribbean regions. 

The GOES-8 image in Figure 25 displays a large region of aerosol dust 
surrounding the Guadaloupe Island region (inside the annotated box). A large cloud mass 
just to the north of Guadaloupe Island region eventually propagated south over the study 
region, which impacted some of the AOD measurements later in the day. 

2. AOD and Phase Function Analysis 

Figures 26 and 27 present the AOD images generated by the NPS algorithm on 
GOES-8 data. The time range for this study is from 1645 UTC through 2045 UTC. The 
data from NOAA-16 at 17:50 UTC and NOAA-14 at 21:16 UTC provided the aerosol 
model index (AMI) to the GOES AOD calculations; these AOD images are displayed in 
Eigure 28. The locations of AOD measurements are shown by the red boxes annotated 
within Eigures 11 and 12. A large cloud field is shown to be propagating southward 
toward Guadaloupe Island throughout the time period. As a result, AOD measurements 
between 19:15 and 20:15 UTC were omitted due to cloud contamination. 

Eigure 29 shows the frequency histograms used to obtain AOD values derived 
from GOES-8 and NOAA-16. AOD values range from 0.33 to 0.48. As the cloud field 
began impacting the Guadaloupe Island region, measurements of AOD were becoming 
increasingly difficult after 18:45 UTC. The normalized distribution pattern shown in 
Eigure 29 became flatter, resulting in larger variability. 

^ NAAPS optical depth values do not directly correspond to the AOD values obtained either by the 
NPS algorithm or AERONET observations. 
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Figure 30 presents the time series of AOD for this case. Both the GOES- 
generated AOD and the AERONET observations of AOD are in good agreement, with 
high AOD values throughout the time period. As mentioned earlier, cloud contamination 
resulted in limited AOD measurements after 19:15 UTC. The lengths of the variability 
bars associated with GOES-8 data increased with time. As a result, the proper selection 
of AOD measurements became increasingly more difficult. The AOD generated from the 
NOAA-16 data at 17:50 UTC (AOD ~ 0.45) and the NOAA-14 data at 21:17 UTC (AOD 
~ 0.48) are also in agreement with AERONET observations. Based on the scatter angle 
profile, local noon occurred toward the beginning of the time period (-16:45 UTC). The 
AOD from AERONET observations tend toward higher values than the satellite-derived 
AOD during the early afternoon hours. The reverse occurs later in the day. 

Eigure 31 displays the phase function value profiles from satellite and AERONET 
data. As in the previous case (Eigure 22), the AERONET phase function values are 
consistently lower than the satellite-derived phase functions. Due to the limited dataset, 
there was not a distinct phase function pattern to describe this dust environment. The 
next chapter will address this issue more clearly. 

C. RESULTS FROM 22 CASES 

Table 1 in Appendix A lists the dates, times and locations of all of the cases in 
this study. The data was selected based on the availability of both satellite and 
AERONET data. Also, entire case study days were rejected when the NOAA data was 
contaminated by sun glint within the region of study. 

Appendix B presents the AOD charts for all 22 cases. Most of the cases have 
AOD profiles that are fairly dust-free (da < 2.5). The exceptions are cases 2, 7, 8, 12, 14, 
16, and 20. Unfortunately, for the high dust cases, data for the analysis is usually limited 
because there tends to be more cloud contamination surrounding all 3 islands within the 
study. As a result the contaminated data are filtered out of the case studies. Within the 
GOES-8 data, local noon occurs during the peak of the scatter angle, between 16:45 and 
17:15 UTC in September. During the late afternoon hours, the AOD profile becomes 
more questionable, as shadows from nearby clouds might affect the clear regions. 
Special care was taken to avoid these problem areas. 
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For each case, measurements were taken in such a way as to compare AOD 
between AERONET and satellite images as close to each other as possible without land 
or cloud contamination. Eigure 32 presents the measurement distances between the 
satellite-derived AOD and the AERONET observations. The easiest region for 
measurements was Bermuda, with an average distance of 19 km between AERONET site 
and the satellite-derived AOD measurement. The islands of Guadaloupe and Puerto Rico 
contain topography that often produces orographic ally enhanced cloud cover on a diurnal 
scale. To avoid cloud contamination, there were several instances where measurements 
had to be deviated by as much as 80 km from the AERONET site. Measurements were 
still conducted in regions further away, so long as the representative AOD value could be 
determined. 

1. Evaluation of the NPS Algorithm 

Eigure 33 compares satellite-derived AOD data with AERONET ‘ground truth’ 
observations for all 22 cases. Within low AOD conditions (5a < 0.2), there is a slight bias 
for NOAA-derived AOD values toward higher values. The GOES-derived values have a 
slight bias toward lower AOD values. Within dust conditions (5a > 0.25), there is a bias 
within the NPS algorithm to an underestimate AOD. Similar findings found by Smith 
(1998) attributed the probable cause to the “no absorption” assumption (tOo = 1) within 
the NPS algorithm. As shown in Eigure 33, the overall standard error for AOD 
measurement is 0.066. Eor the regression analysis, the GOES and NOAA combined 
results in an of 0.67. Individually, the NOAA regression is 0.62 while the GOES 
regression is higher at 0.67. 

Another method to evaluate the performance of the NPS algorithm was to 
determine whether there was bias in AOD calculations due to varying geometries 
between the sun and satellite positions. Eigure 34 is a display that categorizes several 
regression parameters into scatter angle categories. As mentioned earlier, due to the 
position of the study areas and the fixed location of the GOES-8 satellite, local noon 
occurs around the peaks of scatter angles (total backscatter). Therefore, in Eigure 34, the 
scatter angle category of 170“ - 180° (complete backscatter) is within the region of local 
noon. R describes the degree of correlation between satellite-derived AOD and 
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AERONET data. Eor example, = 0.40, indicates that 40% of the original variability of 
the satellite-derived AOD can be explained, with a remaining 60% of residual variability. 
As shown within the bar patterns of Eigure 34, F? values are highest about 140“ - 150“ 
(R^ ~ 0.72) and at 170“ - 180“ (R^ ~ 0.76). Corresponding standard error (S.E.) values are 
at their lowest within the scatter angle categories of 130“ - 140“ and 170“ - 180“, 
respectively. A possible explanation for the higher accuracy about 140 degrees could be 
that the model phase function table of values used within the NFS algorithm converge 
toward one value at -140 degrees. Therefore, there are no aerosol size distribution 
selection errors at this scatter angle. 

Eigure 35 displays the phase function analysis for the 22 cases within dust 
conditions = 0.25). Satellite-derived phase functions (blue dots) and the phase 
functions required to match the AERONET AOD (red dots) are shown. Between 140“ 
and 180“, phase function values generated by the NFS algorithm indicate a pattern of 
higher curvature than that of AERONET-based phase functions. This result is consistent 
with work conducted by Collins et al. (2000) in the ACE-2 experiment off the west 
African coast (upstream from the FRIDE region). Eigure 36 presents their findings. 
Beyond the scatter angle of 140°, non-spherical dust particles were observed to produce a 
flatter phase function shape than non-dust conditions. Eigure 37 shows the 7 aerosol size 
distribution models used in the NFS algorithm with the normali z ed phase function 
patterns supplied by Collins et al. (2000). As shown, within the forward scatter angles 
(0° - 90“), both the spherical and non-spherical contours follow the phase function 
patterns of the 7 models. However, within the backscatter region (90“ - 180“), the non- 
spherical phase function curve deviates from all other curves by revealing a flatter 
profile. As mentioned earlier, the NFS algorithm applies its theory based on non-dust, 
spherical aerosol particles. Based on Figure 37, it would be plausible to apply to the NFS 
algorithm a modified phase function that is flatter in the backscatter region during dust 
events. 
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Figure 11. Composite of SeaWiFs images eentered at 12 UTC on (a) 24 September, 2001 
and (b) 25 September 2001 covering the Atlantic Basin. The location of interest for this 
study is the island of Bermuda. Over the ocean, clear (low aerosol content) regions are in 
dark blue, cloudy regions are solid white, and gray regions depict higher aerosol (dust) 
content. A large plume of dust is visible off of the west coast of Africa. (Courtesy of Dr. 

Douglas L. Westphal at NRL) 
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Figure 12. Plot of NAAPS display of optical depth for sulfate (red shades), dust and smoke 
(green and yellow shades) over the Atlantic Ocean basin for 25 September 2001, 18:00 
UTC. The bottom color bar shows the AOD range for dust and smoke. 850 mb model¬ 
generated wind barbs are also displayed. (Courtesy of Dr. Douglas L. Westphal at NRL) 
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Figure 13. GOES-8 visible image on 25 September 2001 at 17:15 UTC. The annotated box 

surrounds the region of Bermuda. 
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Figure 14. Time series of AOD images generated for 25 September 2001 from GOES-8 data 
that surrounds Bermuda. Pixel sizes are 1.1 km by 1.1 km and the domain is 
approximately 110 km by 110 km. The times range from 13:15 UTC to 16:15 UTC. Red 
boxes depiet loeations where the representative AOD for that area was measured. AOD 
eolor eontours are defined on the left side of eaeh image. 
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Figure 15. Continuation of Figure 14 for times ranging from 17:15 UTC to 20:15 UTC. 
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Figure 16. AOD image generated for 25 September 2001 at 18:17 UTC from NOAA-16 data 
that surrounds Bermuda. Pixel sizes are 1.1 km by 1.1 km and the domain is 
approximately 110 km by 110 km. Red box depicts the location where the representative 
AOD for that area was measured. AOD color contours are defined on the left side of the 

image. 
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Figure 17. An example of a histogram frequency plot for determining AOD. A normalized 
curve (solid red curve) displays the distribution, where the center of mass (*) is situated 
at the top of the curve. The resulting AOD value of 0.03 is then obtained. The variability 
measurements of AOD are obtained from a line drawn at half the amplitude of the 
normalized curve that indicates one standard deviation from the measured AOD value. 
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Figure 18. Time series of histogram plots of AOD for 25 September 2001, between 13:15 
UTC and 16:45 UTC. Plots are generated from boxes defined in Figure 14. The AOD 
values are shown in red, while the two adjacent values represent the variability of AOD 
measurements. The red annotations depict the AOD measurement process that is 

described in Figure 17. 
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Figure 19. Continuation of time series of histogram plots from Figure 15 of AOD for 25 
September 2001, between 17:15 UTC and 20:15 UTC. 
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Optical depth 

Figure 20. Histogram plot of AOD for 25 September 2001 at 18:17 UTC from NOAA-16 
data. The plot is generated from box defined in Figure 16. Times are in UTC. The AOD 
values are shown in red, while the two adjacent values represent the variability in the 

AOD measurement. 
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AERONET observations (yellow) for 25 September 2001 over Bermuda. NOAA data 
are presented as individual blue dots whereas GOES-8 data are connected with a blue 
line. Red bars represent variability in AOD measurements. The corresponding scatter 
angle profile is represented as a purple dashed line. 
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Scatter angle (degrees) 

Figure 22. Plot of phase function as a function of scatter angle for the case of 25 September 
2001 over Bermuda. Blue diamonds indicate phase functions obtained from satellite 
data. Red boxes indicate phase function values required for the satellite-derived AOD to 

match AERONET’s AOD. 
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Figure 23. Composite of SeaWiFs images centered at 12 UTC on (a) 17 September 2001 and 
(b) 18 September, 2001 covering the Atlantic Basin. The location of interest for this 
study is Guadaloupe Island. Over the ocean, clear (low aerosol content) regions are in 
dark blue, cloudy regions are solid white, and gray regions depict higher aerosol (dust) 
content. A large plume of dust is visible off of the west coast of Africa. (Courtesy of Dr. 

Douglas L. Westphal at NRL) 
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Figure 24. Plot of NAAPS display of optical depth for sulfate (red shades), dust and smoke 
(green and yellow shades) over the Atlantic Ocean basin for 18 September 2001, 18:00 
UTC. (Courtesy of Dr. Douglas L. Westphal at NRL) 


49 



















































Figure 25. GOES-8 visible image on 18 September 2000 at 17:15 UTC. The annotated box 

surrounds the region of Guadaloupe Island. The area of aerosol dust is also annotted. 
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Figure 26. Time series of AOD images generated for 18 September 2001 from GOES-8 data 
that surrounds Guadaloupe Island. Pixel sizes are 1.1 km by 1.1 km and the domain is 
approximately 110 km by 110 km. The times range from 16:45 UTC to 18:15 UTC. Red 
boxes depiet locations where the representative AOD for that area was measured. AOD 
color contours are defined on the left side of each image. 
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Figure 27. Continuation of Figure 26 for times ranging from 19:15 UTC to 20:45 UTC. 
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Figure 28. AOD images generated for 18 September 2001 at 17:50 UTC and 21:16 UTC 
from NOAA-16 and NOAA-14 data, respectively. Images surround Guadaloupe Island 
Pixel sizes are 1.1 km by 1.1 km and the domain is approximately 110 km by 110 km. 
Red boxes depict the locations where the representative AOD for that area were 
measured. AOD color contours are defined on the left side of each image. 
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Figure 29. Time series of frequency histogram plots of AOD for 18 September 2001, 
between 16:45 UTC and 21:16 UTC. The top 6 plots are derived from GOES-8 data and 
the bottom 2 plots are derived from NOAA-16 and NOAA-14 data, respectively. Plots 
are generated from boxes defined in Figures 26 and 27. For each plot, the AOD value is 
shown in red, while the pair of black values represent the variability of AOD 
measurements. The red lines and curves depict the AOD measurement process that is 

described in Figure 17. 
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Figure 30. Chart of AOD for both satellite-derived data (navy blue dots) and AERONET 
observations (yellow) for 18 September 2001 over Guadaloupe Island. NOAA data are 
presented as individual blue dots whereas GOES-8 data are connected with a blue line. 
Red bars represent variability in AOD measurements. The corresponding scatter angle 
profile is represented as a purple dashed line. 
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Scatter angle (degrees) 


Figure 31. Plot of phase function as a function of scatter angle for the case of 18 September 
2001 over Guadaloupe Island. Blue diamonds indicate phase functions obtained from 
satellite data. Red boxes indicate the phase function values required for the satellite- 
derived AOD to match AERONET AOD. 


56 




















a)»■ 

u 

S <0 ■ 

P 

IS) M - 

■H 

Q „ 


Bermuda 


r. . 






Avg dist: 19 km 
74 Observations 


Guadaloupe Island 




Figure 32. Distribution of distance measurements between each of the 3 AERONET sites 
and the location of the satellite-derived AOD for all 22 cases. Eor each AERONET 
region, the average distance from the AERONET station to the center of the satellite- 
derived measurement box is annotated in red. 
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AERONET-derived AOD 


Figure 33. Comparisons between NPS-calculated AOD and AERONET-observations of 
AOD for the 22 cases as listed in Appendix A.l. The red line depicts one to one 

relationship. 
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Figure 34. Evaluation of NFS algorithm, partitioned over eategories of seatter angles. 
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Figure 35. Comparisons of phase funetions between satellite data and AERONET data for 22 
cases. Eigure only depicts dust conditions (AERONET AOD = 0.25). 



Eigure 36. Phase function plots for the free troposphere, dust layer obtained from the second 
Aerosol Characterization Experiment (ACE-2). Data is supplied by Collins, et al. (2000). 

Phase function values are normalized. 
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Figure 37. Phase function plots comparing the dust plots in Collins et al. (2000) (thick 
curves) with the model phase functions used in the NPS algorithm (thin curves). 
Annotated numbers are used to identify the model numbers. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

A radiative transfer algorithm (NPS algorithm) that processes AOD within a 
cloud-free maritime atmosphere was presented for the NOAA POES AVHRR and GOES 
imager. This algorithm applies linearized, single-scatter theory with an estimate of bi¬ 
directional surface reflectance. Using a technique developed by Durkee et al. (1991), 
scattering phase functions are parameterized to seven aerosol size distributions by 
applying the ratios of channel 1 and 2 radiances (S 12 ) of the AVHRR. The S 12 value is 
then translated to an aerosol model index (AMI) value that is accessed by GOES 
processing. Both the NOAA and GOES processing then apply the AMI to a look up table 
to determine the scattering phase function. The development and validation of the NPS 
algorithm is a continuation of the initial work performed by Brown (1997). 

Unlike the AVHRR data, the GOES-8 visible radiance data required an 
unconventional calibration scheme developed by Rao et al. (1999) and Rao and Zhang 
(1999). During the analysis period of this paper (July and August, 2001) signal gain 
factors ranging from 1.727 during the first case study (July 27, 2001) through 1.740 
during the last study (September 25, 2001) were applied to the GOES raw channel 1 
albedo data. In addition, preliminary comparisons of AOD between GOES and NOAA, 
and AERONET revealed that an additional GOES-8 channel 1 correction factor of 
-5.5Wm'^ sr‘was necessary. 

To validate the NPS algorithm, retrieved AOD data was collected from 22 cases 
during July and September of 2001 and compared to AERONET radiometers positioned 
within 3 island locations of the western Atlantic Ocean. Eor each case, a time series 
format was used to study temporal variations in AOD. Overall, the comparison dataset 
has a correlation coefficient of 0.67 with a standard error of 0.07. Within higher AOD 
cases (da > 2.5), the general trend was for the NPS-generated AOD values to under¬ 
estimate AERONET-observed conditions, probably due to assumptions of non-absorption 
applied to aerosol particles. When related to scatter angles, the NPS-generated AOD 
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calculations performed best within the backscatter angle ranges of 130“ - 140° and 170° - 
180°. 


A major part of the uncertainty to the AOD processing is the proper selection of 
the scattering phase function. The problem with the AOD results generated from GOES 
data is that the aerosol conditions over a particular region were assumed to be unchanging 
throughout the entire time period, which might last for 8 hours. Unfortunately, trying to 
validate a representative aerosol distribution was beyond the scope of this study. The 
AOD results were apparently not very sensitive to any aerosol model assigned to the 
calculation. 

The pattern of the NFS algorithm- generated phase functions was evaluated during 
the dust events (da = 0.25). Between the scatter angles of 140° through 180°, the phase 
function pattern required to match AERONET observations indicated a flatter profile 
than that produced by the NFS algorithm on GOES and NOAA data. The AERONET- 
generated phase function profile was confirmed with studies by Collins et al. (2000) for 
non-spherical dust aerosols compared to spherical particles assumed in the NFS 
algorithm. 

B. RECOMMENDATIONS 

Several research initiatives would aid in the improvement of the NFS algorithm: 

• Coordinate with aerosol modelers to develop a more accurate assessment 
of aerosol conditions on a global scale. The NAAFS model developed at 
NRE is one such source. The NAAFs would greatly help the NFS 
algorithm in determining the proper phase function parameters. 

• Apply the NFS algorithm to other satellite platforms. Freliminary studies 
by Brown (2001) indicate that the SeaWiEs sensor, with its suite of visible 
sensing channels, could greatly enhance the S 12 phase function 
parameterization technique. Another instrument with great potential 
applications is the Moderate Resolution Imaging Spectroradiometer 
(MODIS). 

• Improve sun glint removal process by incorporating wind speed, which is 
available from the Defense Meteorological Satellite Frogram’s (DMSF) 
Special Sensor Microwave/Imager (SSM/I) data as well as from a variety 
of regional and global numerical models such as COAMFS^'^ and 
NOGAFS. 
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Develop an on-board calibration to the GOES channel 1 sensors that will 
greatly reduce AOD uncertainties. 

Refine the phase functbn lookup table to include more realistic aerosol 
size distributions for dust events. 

Apply appropriate single-scatter albedo (? o) values for conditions with 
known absorption properties. 
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APPENDIX A 


Table 5. List of 22 cases that include the GOES and NOAA pass times (in UTC) and 
locations. The case studies are identified in sequential numbers, shown in parentheses. 


Dates 

Satellite Times 

Study iocations 

GOES corr. factors 

(1) 

GOES-8: 16:45-21:15 



27 July 

NOAA-16: 18:42, NOAA-14: 21:51 

NOAA-14: 21:51 

Guadaloupe Island 

1.727 

(2) 

GOES-8: 17:15-21:15 



28 July 

NOAA-16:18:31,NOAA-14 :21:39 

La Paguera 

1.727 

(3) 

GOES-8: 17:15-21:15 



02 Sept 

NOAA-16: 17:18, NOAA-14: 21:08 

Bermuda 

1.734 

(4) 

GOES-8: 16:45-21:15 



04 Sept 

NOAA-16: 18:35 

Guadaloupe Island 

1.735 

(5) 

GOES-8: 17:15-21:15 



05 Sept 

NOAA-16: 18:26 

La Paguera 

1.735 

(6) 

GOES-8: 17:15-21:15 



06 Sept 

NOAA-16: 18:16 

Guadaloupe Island 

1.735 

(7) 

GOES-8: 17:15-21:15 



07 Sept 

NOAA-16: 18:05 

Guadaloupe Island 

1.735 

(8) 

GOES-8: 17:15-21:15 



08 Sept 

NOAA-16: 17:55 

Guadaloupe Island 

1.736 

(9) 

GOES-8: 15:15-20:45 



11 Sept 

NOAA-14: 21:00 

Bermuda 

1.736 

(10) 

GOES-8: 15:15-20:45 



11 Sept 

NOAA-14: 21:00 

Guadaloupe Island 

1.736 

(11) 

GOES-8: 16:15-20:45 




17 Sept 


NOAA-16: 18:01 


Bermuda 


1.737 

















































(12) 

GOES-8: 16:15-20:45 



17 Sept 

NOAA-16: 18:01 

Guadaloupe Island 

1.737 

(13) 

GOES-8: 16:45-21:15 



18 Sept 

NOAA-16: 17:50, NOAA-14: 21:16 

Bermuda 

1.738 

(14) 

GOES-8: 16:45 - 20:45 



18 Sept 

NOAA-16:17:50, NOAA-14: 21 :16 

Guadaloupe Island 

1.738 

(15) 

GOES-8: 15:15-21:15 



19 Sept 

NOAA-16: 17:40, NOAA-14: 21:04 

Bermuda 

1.738 

(16) 

GOES-8: 15:15-18:45 

La Paguera 


19 Sept 

NOAA-16: 17:40 


1.738 

(17) 

GOES-8: 14:45-21:15 

Bermuda 


20 Sept 

NOAA-16: 17:29, NOAA-14: 20:53 


1.738 

(18) 

GOES-8: 14:15-21:15 



21 Sept 

NOAA-16: 17:19, NOAA-16: 19:00 

La Paguera 

1.738 

(19) 

GOES-8: 14:15-21:15 



22 Sept 

NOAA-16: 18:47 

La Paguera 

1.738 

(20) 

GOES-8: 17:15-21:15 



23 Sept 

NOAA-16: 18:38 

La Paguera 

1.739 

(21) 

GOES-8: 13:15-20:45 



25 Sept 

NOAA-16: 18:17 

Bermuda 

1.739 

(22) 

GOES-8: 13:15-21:15 



29 Sept 

NOAA-16: 17:35 

La Paguera 

1.740 
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APPENDIX B. 


Time series of AOD for 22 cases ranging from 27 July, 2001 through 29 
Spetember, 2001 over the regions of Bermuda, Guadaloupe Island, and La Paguera, 
Puerto Rico. Both satellite-derived data (navy blue dots) and AERONET observations 
(yellow) NOAA data are presented as individual blue dots whereas GOES-8 data are 
connected with a blue line. Red bars represent wiability in AOD measurements. The 
corresponding scatter angle profile is represented as a purple dashed line. 


Case No. 1 20010727 Region: Guadaloupe Island 
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Case No. 2 20010728 Region: La Paguera, Puerto Rico 



Case No. 3 20010902 Region: Bermuda 
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Case No. 4 20010904 Region: Guadaloupe Island 



Case No. 5 20010905 Region: La Paguera, Puerto Rico 
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Case No. 6 20010906 Region: Guadaloupe Island 



Case No. 7 20010907 Region: Guadaloupe Island 



72 
















































Case No. 8 20010908 Region: Guadaloupe Island 



Case No. 9 20010911 Region: Bermuda 



Time (UTC) 
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Case No. 10 20010911 Region: Guadaloupe Island 



Time (UTC) 


Case No. 11 20010917 Region: Bermuda 
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Case No. 12 20010917 Region: Guadaloupe Island 



Case No. 13 20010918 Region: Guadaloupe Island 
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Case No. 14 20010918 Region: Guadaloupe Island 



Case No. 15 20010919 Region: Bermuda 
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Case No. 16 20010919 Region: La Paguera, Puerto Rico 



Case No. 17 20010920 Region: Bermuda 
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Case No. 18 20010921 Region: La Paguera, Puerto Rico 



Case No. 19 20010922 Region: La Paguera, Puerto Rico 
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Case No. 20 20010923 Region: La Paguera, Puerto Rico 



Case No. 21 20010925 Region: Bermuda 
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Case No. 22 20010929 Region: La Paguera, Puerto Rico 
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